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^ ■ A superposition rule is a particular type of map that enables one to express the 

general solution of certain systems of first-order ordinary differential equations, the so- 
CO ' called Lie systems, out of generic families of particular solutions and a set of constants. 

The first aim of this work is to propose several generalisations of this notion to second- 
order differential equations. Next, several results on the existence of such generalisations 

r^ I arc given and relations with the theories of Lie systems and quasi-Lic schemes are found. 

C^H' Finally, our methods are used to study second-order Riccati equations and other second- 

(-H \ order differential equations of mathematical and physical interest. 
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The origin of the study of superposition rules can be traced back to 1893, when Lie, Guld- 
berg, and Vessiot published a series of papers [28, 34, 41, 42] characterising and analysing 
^^ ' systems of first-order ordinary differential equations admitting a superposition rule. This 

^P ■ theory, broadly analysed during its beginning, was rarely treated for the next eighty years. 

Nevertheless, the interest in this topic has revived in the last few decades, and many works 
have recently been devoted to the analysis of its properties and to the study of its applica- 
/\^ • tions and generalisations [2, 4, 6, 10, 12, 13, 25, 29, 32, 33, 39, 43, 45]. Among these works, 

we can emphasise the results described in [6, 29, 45]. 

In spite of its interesting properties and applications in Mathematics, Control Theory, 
and Physics (see [6, 9, 10, 12, 13, 25, 36, 39]), the theory of Lie systems possesses a certain 
lack of applicability in many other problems described by differential equations that are not 
Lie systems. This motivates the generalisation of the methods for analysing Lie systems 
carried out here, and in previous works, so as to investigate non-Lie systems. 

The theories of Lie systems [6, 34, 45] and quasi-Lie schemes [4, 8] mainly concern the 
investigation of systems of first-order differential equations. Nevertheless, various attempts 
have been carried out to apply their methods to study second-order differential equations 
(SODEs) [9, 10, 12, 25, 39, 43]. These papers were rather more focused on applying the 
methods of the aforementioned theories to SODEs than on describing new theoretical results. 
This explains why, although these applications suggested the existence of a new type of 
(nonlinear) expressions describing general solutions of SODEs, no further discussion about 
such expressions and their properties was carried out. 



In view of the above comments, the first aim of this work is to define some new notions 
describing, as particular cases, those expressions for the study of SODEs appearing in the 
very recent hterature, namely, time-dependent and time-independent superposition rules for 
systems of SODEs. Furthermore, some theoretical results concerning the existence of such 
superposition rules are demonstrated, and several relations with the theories of Lie systems 
and quasi-Lie schemes are shown. In this way, our achievements allow us to clarify diverse 
procedures and techniques used in previous works to investigate SODEs. 

After developing the theoretical part of our work, our results are applied to study families 
of SODEs appearing broadly in the Physics and Mathematics literature [1, 7, 14, 15, 18, 
25, 30, 31, 38]. More specifically, time-independent and time-dependent superposition rules 
are derived for such families using the theory of Lie systems and quasi-Lie schemes. This 
reduces the determination of the general solution of any instance of that families to the 
determination of four of its particular solutions. On one hand, this provides a modern fully 
geometric demonstration for a result sketched in a Vessiot's work [43] that was derived in 
ad hoc way. On the other hand, our work supplies a method to analyse the solutions of 
specific forms of time-dependent Lienard equations [15, 25, 38], various interesting equations 
studied by Ince and Davis [18, 30], some modified Emden equations [38], certain second- 
order Riccati equations [1, 7, 14], etc. Moreover, particular forms of the analysed SODEs 
arise in the analysis of fusion pellets [21]. 

Among the aforementioned applications, particular attention is paid to the study of 
second-order Riccati equations. The interest is mainly due to two reasons. On one hand, 
these equations are widely studied in order to analyse their mathematical properties [1, 14, 
17, 20, 23, 27, 35]. On the other hand, second-order Riccati equations and some of their 
particular cases, e.g. certain modified Emden equations [15, 16, 38] or various Painleve-Ince 
equations [3, 14, 21, 26, 31], appear in many physical problems, like the problem of waves 
on shallow water. 

Another application of second-order Riccati equations deserves a special attention: these 
equations are an element of the so-called Riccati chain describing the Backlund transforma- 
tions for various partial differential equations [7, 27]. More specifically, for each one of such 
PDEs, a particular solution defines a one-parametric family of n-order Riccati equations. 
The solutions of the members of such a family allow one to build up a new family of solutions 
of the initial PDE [27]. Obviously, the time-dependent superposition rule derived in this 
work for second-order Riccati equations simplifies the study of the solutions of such a family 
and, in consequence, the determination of new solutions for those PDEs whose Backlund 
transformations are determined by these equations. 

The plan of the paper is as follows. Section 2 mainly recalls the notions of the theory 
of Lie and quasi-Lie systems necessary to follow the results of our paper. Additionally, the 
usefulness of quasi-Lie schemes to explicitly determine time-dependent superposition rules 
for systems of first-order differential equations is discussed. Section 3 is concerned with one 
of the novelties of our paper: the definition and analysis of the time-independent and the 
time-dependent superposition rule notions for systems of SODEs. In Section 4 our previous 
theoretical results are employed to derive a superposition rule for all the elements of a 
family of SODEs with many applications to Physics and Mathematics. A quasi-Lie scheme 
is used in Section 5 to derive a time-dependent superposition rule for second-order Riccati 
equations. Finally, some calculations to prove various results of the paper are detailed in 
the Appendix. 



2 Lie systems and quasi-Lie schemes 

Here we mainly review some features about Lie systems, quasi-Lie schemes and time- 
dependent superposition rules for systems of first-order differential equations [4, 5, 6, 8]. 
In addition, the use of quasi-Lie schemes for deriving time-dependent superposition rule is 
discussed. For simplicity, we restrict ourselves to analysing systems of differential equations 
on linear spaces and we skip various technical details that are not relevant to understand 
our procedures. 

Lie initiated in [34] the nowadays called theory of Lie systems when started investigating 
the conditions that ensure that a system of first-order ordinary differential equations of the 
form 

^=X\t,x), i = l,...,n, (1) 

admits a superposition rule, i.e. a t-independent function $ : R"™ x M" — )• R", x = 
$(xi, . . . , Xm', Ai, . . . , A„), such that the general solution of the system can be written as 

X{t) = ^[xi{t),...,Xm{t);\l,...An), (2) 

where {xa{t) \ a = 1, . . . , m} is any generic family of particular solutions and Ai, . . . , A^ is 
a set of n constants. 

Note that the general solution of a linear homogeneous system of first-order differential 
equations in M cannot be cast into the form (2) for every set of n particular solutions: they 
must be linearly independent. In a similar way, for all known Lie systems (cf. [32, 45]), 
their expressions (2) only hold for certain families of particular solutions. More specifically, 
it is said that expression (2) is valid for any 'generic' family of m particular solutions if 
there exists an open dense subset U C M"'" such that expression (2) is satisfied for every 
set of particular solutions xi(t), . . . , Xm{t) such that (xi(0), . . . , Xm{^)) lies in U . Obviously, 
'almost every' set of m particular solutions holds the above property and this is, indeed, the 
approximate meaning of the term 'generic' in the above definition. 

Lie found a characterisation describing systems of first-order differential equations admit- 
ting a superposition rule [34] and such a characterisation was recently reformulated in terms 
of the modern language of Differential Geometry [6]. In these terms, each system of the form 
(1) is described by means of a time-dependent vector field X{t, x) = Yll=i X^i't, x)d/dx^ on 
M". This description is used to establish that a first-order system (1) admits a superposi- 
tion rule if and only if its associated time-dependent vector field X(t, x) can be written as 
a linear combination 

r 

X{t,x) = Y,bait)X^ix), (3) 



a=l 



where the vector fields 



" d 



are a basis for an r-dimensional real Lie algebra V of vector fields, the so-called associated 
Vessiot-Guldberg Lie algebra. In other words, a time-dependent vector field X{t, x) describes 
a Lie system if and only if the family of vector fields {Xt}t^R, with Xf : x £ M" — )• Xt{x) = 



X{t, x) G TM", holds {Xt}t£M. C V^ for a certain finite-dimension Lie algebra of vector fields 
V. 

Various procedures are available to derive a superposition rule. Let us sketch here one of 
these methods to be used within this work (for a full description, see [6]). The key element 
of this procedure is the so-called diagonal prolongation of a vector field. Given a vector field 
X{x) = ^^=1 X^ {x)d / dx'^ over M", its diagonal prolongation to (IR"')^"*"^ is the vector field 
over this space 

p n a 



X{xo,... ,Xm) = 7 ./^^'(^Q 



a=0 i=l ^^"^ 

Now, our approach to derive a superposition rule for system (1) starts by determining a 
decomposition (3) for its associated time-dependent vector field. Next, it must be deter- 
mined the natural number m so that the diagonal prolongations to (M")™ of the vector fields 
Xi, . . . ^Xj. become linearly independent at each point of an open dense subset of this space. 
Note that the diagonal prolongations Xi, . . . , Xr to M'^(™-+^) of the vector fields Xi, . . . , Xj- 
satisfy their same commutation relations and they are linearly independent over an open 
dense set U C M'^('"+^) again. Consequently, they span an involutive distribution 

Vs = {Xi[x),...,Xr{x)), X£U, 

of rank r over U. The vector fields of this distribution admit n{m + 1) — r common local 
first-integrals. Among them, one can choose n first-integrals giving rise to a n-codimensional 
local foliation T horizontal with respect to the projection n : {xq, . . . ,Xm) £ (M")™+^ i— )• 
(xi, . . . , Xm) S (M")'". Roughly speaking, the former first-integrals enable us to express the 
coordinates Xq, with i = 1, . . . ,n, in terms of the other variables and the n first-integrals, 
giving rise to the superposition rule for system (1). 

Since it would be interesting to find a way to apply the methods for analysing Lie 
systems to a broader set of differential equations, the theory of quasi-Lie schemes, whose 
basic features are described below, was developed [4]. Although such a theory applies 
to systems of first-order differential equations associated with complete and non-complete 
time-dependent vector fields, the forthcoming presentation will focus on studying systems 
associated with complete vector fields in order to simplify the presentation and underline 
the main results of the theory. Nevertheless, a careful analysis shows that our main claims 
remain valid for systems associated with non-complete vector fields, although with some 
technical minor modifications. 

Let us now turn to define the key notion of the theory of quasi-Lie schemes. 

Definition 2.1. A quasi-Lie scheme S{W,V) is made up by two finite-dimensional vector 
spaces of vector fields W, V satisfying the following conditions: 

• W is a linear subspace of V. 

• W is a Lie algebra of vector fields, that is, [W, W] C W. 

• W normahses V, i.e. [W, V] CV. 

Roughly speaking, each quasi-Lie scheme S{W, V) determines two families of time- 
dependent vector fields, V(M) and VF(M), which are made of time-dependent vector fields 
taking values in V and W, respectively. The first family, 1^(M), is intended to describe the 



first-order systems which can be analysed by means of S{W,V). The second one, VF(M) 
is used to define a group of time-dependent changes of variables: the group of the scheme 
G{W). Now, each system related to a time-dependent vector field of V{M) transforms, under 
any element of Q(W), into another system related to another element of y(]R). This is the 
most important property of quasi-Lie schemes [4, Proposition 1]. It has been widely used in 
recent works [4, 8, 11] to transform diverse types of differential equations, e.g. Abel equa- 
tions [11], into other differential equations of the same type and, more specifically, into Lie 
systems. In this way, the theory of Lie systems applies to investigate this last system and, 
by undoing the performed change of variables, multiple properties of the original system 
can be stated. 

Let us explain the above claims more carefully. Every time-dependent vector field X 
over M" gives rise to a generalised flow g^ , i.e. a map g^ : {t,x) G M x R" i— )• gf{x) = 
g^{t,x) G R"", with g^ = IdiRn, satisfying that 7^(t) = gf{xo) is the integral curve of the 
time-dependent vector field X passing through the point xq € R" at i = 0. Now, since W 
is a Lie algebra of vector fields, it can be proved that the generalised fiows corresponding 
to time-dependent vector fields with values in W form a group, the so-called group of the 
scheme Q(W), with composition law {g -k h)t = gt o ht and g,h £ G{W). The neutral 
element is e : (t, x) € R x R" i— ;• x G R" and every generalised flow g admits an inverse 
g~^ : {t,x) G R X R" i— )• [gt)~^{x) G R". As every generalised flow can be considered 
as a time-dependent change of variables, the group G{W) can be regarded as a group of 
time-dependent changes of variables. 

Generalised flows can also act on time-dependent vector fields. Given a time-dependent 
vector field Y and a generalised fiow h, the action of h over X, let us say h^X, is the time- 
dependent vector field whose generalised fiow is h-kg . In this terminology, the main result of 
the theory of quasi-Lie schemes [4, Proposition 1] establishes that for every time-dependent 
vector field Y of V(R) and g G G{W), the time-dependent vector field g^Y belongs to F(R). 
In other words, the family of systems associated with the elements of V{M) is stable under 
the time-dependent changes of variables of G(W). Among such systems, special relevance 
have those ones, the so-called Lie systems, that can be transformed into Lie systems. The 
precise definition of this notion is detailed below. 

Definition 2.2. A system of differential equations describing the integral curves of a time- 
dependent vector field X is a quasi-Lie system, with respect to a quasi-Lie scheme S{W, V), 
if there exist a g £ G{W) and a Lie algebra Vq C V such that g-f^X is a time-dependent 
vector field taking values in Vq. 

Quasi-Lie systems admit certain properties as a consequence of its relation to Lie sys- 
tems. For instance, each quasi-Lie system admits its general solution to be described in 
terms of each generic family of particular solutions, a set of constants, and the time, i.e. it 
admits a time- dependent superposition rule [4]. Moreover, the theory of quasi-Lie systems 
provides powerful methods to explicitly determine such superpositions [4, 8, 11]. In order 
to understand the relevance of these facts, it is necessary to pay attention to the following 
remarks. 

Every system of first-order differential equations admits time-dependent superposition 
rules [4, 5]. For instance, every first-order system describes the integral curves of a time- 
dependent vector field admitting a generalised flow which is indeed a particular type of 
time-dependent superposition rule for the system. In spite of this, determining such time- 



dependent superposition rules can be, like in the case of the aforementioned example, as dif- 
ficult as solving the initial system [4, 5]. Thus, what it really matters about time-dependent 
superposition rules is the description of procedures to determine them explicitly as, for 
instance, quasi-Lie schemes. 

Apart from the above remark, there is another reason to use quasi-Lie schemes to de- 
termine time-dependent superposition rules: every quasi-Lie scheme provides a family of 
first-order systems admitting the same time-dependent superposition rule as a bonus [5, 
Proposition 14]. Let us briefly analyse this fact. Given a quasi-Lie system with respect to 
a quasi-Lie scheme S(W,V), there exists an element g G Q{W) and a Lie algebra Vq C V 
such that the time-dependent vector field X associated with the quasi-Lie system holds 
that g-tfX takes values in Vq. Therefore, the set Sg{W,V;Vo) of quasi-Lie systems with 
respect to S{W, V) satisfying that g^X takes values in Vq is not empty. Moreover, it is 
made of time-dependent vector fields of the form X' = {g~^)^Y, with Y being any time- 
dependent vector field taking values in Vq. As all the elements of Sg(W,V;Vo) admit a 
common time-dependent superposition rule (cf. [5, Proposition 14]), we have determined a 
family of systems admitting the same time-dependent superposition rule than X. 

3 Superposition rules for systems of SODEs 

Previous studies about second-order differential equations carried out by means of the theory 
of Lie systems and quasi-Lie schemes were lacking in a theoretical explanation of the methods 
and notions employed [4, 9, 10, 39, 43]. The main aim of this section is to provide a definition 
for the new notions appearing, with no further explanation, within previous works along with 
a theoretical explanation of the methods there performed. In this way, a basic background 
for the posterior theoretical treatment of the subject is laid down. 

Recall that the theory of Lie systems was initiated as a result of the study of systems of 
first-order differential equations admitting their general solutions to be expressed in terms 
of each generic family of particular solutions and a set of constants. Nevertheless, not only 
such systems admit their general solutions to be described in this way. For instance, every 
second-order differential equation of the form x = a(t)x, with a{t) any time-dependent 
function, satisfies that its general solution, x{t), can be put into the form 

x{t) = \iXi{t)+X2X2{t), (4) 

with Ai, A2 being two real constants and xi{t),X2{t) being any family of two particular solu- 
tions such that {xi{t),xi{t)) and {x2{t),X2{t)) are, for every t G M, two linearly independent 
elements of TM ~ M^. In a similar way, other expressions determining the general solution 
of certain SODEs have recently been described in the literature [9, 12]. This suggests us to 
propose the following definition which covers, as particular cases, all the previous expressions 
occurring in the literature. 

Definition 3.1. We say that a system of second-order differential equations 

x" = F\t,x,x), i = l,...,n, (5) 

on M" admits a superposition rule if there exists a map ^ : TM™" x M?"^ — t- M" such that its 
general solution, x{t), can be written as 

x{t) = '^{xl(t),...,Xr,^(t),Xl{t),... , Xmit)] Xl, ■ ■ ■ , X2n), (6) 

6 



in terms of each generic family, xi{t), . . . ,Xmit), of particular solutions, their derivatives, 
and a set of 2n constants. 

In order to grasp the previous definition, it is necessary to precisely establish the meaning 
of 'generic' in the above statement. Formally, it is said that expression (6) is valid for a 
generic family of particular solutions when it holds for every family of particular solutions, 
xi{t), . . . ,Xm{t), satisfying that (xi(0),i;i(0), . . . ,Xm(0),Xm(0)) G U, with U being an open 
dense subset of TM"™. In this way, as in the case of superposition rules for Lie systems, the 
term 'generic' amounts to 'almost every'. 

In order to note that all those aforementioned expressions studying the general solutions 
for certain SODEs are superposition rules in the above sense, it is necessary to explain 
an important detail. Some of such expressions, like the one for studying Milne-Pinney 
equations, depend on a generic set of m particular solutions, a set of constants and a set 
of time-independent constants of the motion [9]. These expressions seem to differ from the 
above definition. Nevertheless, if we take into account that such constants of the motion 
depend on the m particular solutions, we can notice that, indeed, they are superposition 
rules in the above sense. 

Although there exists no characterisation for systems of SODEs of the form (5) admitting 
a superposition rule, there exists a special class of such systems, the so-called SODE Lie 
systems, accepting such a property. Even though this fact has been employed implicitly in 
the literature, it has never been proved explicitly. In view of these facts, we next furnish the 
definition of the SODE Lie system notion along with a proof showing that every SODE Lie 
system admits a superposition rule. In addition, some remarks about the interest of this 
notion and its main properties are discussed. 

Definition 3.2. We say that the system of SODEs (5) is a SODE Lie system if the system 
of first-order differential equations 

{X* = v\ 
v' = F\t,x,v), 

obtained by adding the new variables v^ = i*, with i = 1, . . . ,?7-, to system (5), is a Lie 
system. 

Proposition 1. Every SODE Lie system, (5) adm,its a superposition rule ^ : TM""* xM^" — )• 
W of the form ^ = vr o $, where <& : TM"'" x M^** — )• TM" is a superposition rule for the 
system (1) and vr : TM" — )• M" is the projection associated with the tangent bundle TM". 

Proof. Each SODE Lie system of the form (5) is associated with a first-order system of 
differential equations (7) admitting a superposition rule $ : TM"™ x R^" — >• TM"-. This 
allows us to describe the general solution (x(t), v{t)) of system (7) in terms of each generic 
set {xa{t),Va{t)), with a = 1, . . . , 7TT., of its particular solutions and a set of 2n constants, i.e. 

{x{t),v{t)) = $ (Xi(t), . . . , X^{t),Vi{t), ..., Vm{t);Xl, . . . , A2n)- (8) 

Obviously, each solution, Xp{t), of the second-order system (5) corresponds to one and only 
one solution {xp{t),Vp{t)) of (1) and viceversa. Since {xp{t) , Vp{t)) = {xp{t),Xp{t)), it follows 
that the general solution x{t) of (5) can be put in the form 

x{t) = vr O $ {Xi(t), ..., Xm(t),Xi{t), ..., Xm{t)] Ai , . . . , \2n) , (9) 



where Xa{t), with a = 1, . . . ,n, is a generic family of particular solutions of (5). In other 
words, the map ^ = vr o $ is a superposition rule for the system of second-order differential 
equations (5). D 

Since every autonomous system is related to a one-dimensional Vessiot-Guldberg Lie 
algebra [4], it straightforwardly follows, from the above proposition, next corollary. 

Corollary 1. Every autonomous system of second- order differential equations of the form 
X* = F^{x, x), with i = 1, . . . ,n, admits a superposition rule. 

Despite its theoretical interest, the above result cannot straightforwardly be used to 
derive superposition rules most of the times. Actually, the superposition rule guaranteed by 
Proposition 1 relies on obtaining a superposition rule for an autonomous first-order system 
of differential equations. Considering the method sketched during Section 2, we find that 
determining this superposition rule demands the determination of all the integral curves of 
a vector field on (TM*^)^. Although the existence of the solution of this problem is known, 
its explicit description can be as difficult as solving the initial system (indeed, this is usually 
the case). Consequently, deriving explicitly a superposition rule for the above autonomous 
system frequently relies on searching an alternative superposition rule for the associated 
first-order system. 

Many superposition rules for systems of second-order differential equations do not show 
an explicit dependence on the derivatives of the particular solutions. Consider, for instance, 
either the linear superposition rule (4) for the equation x = a{t)x, or the affine one, 

X{t) = Xl{xi{t) - X2{t)) + \2{x2{t) - Xs{t)) + X3{t), 

for X = a{t)x + h{t). Such superposition rules are called velocity free superposition rules or 
even free superposition rules. The conditions ensuring the existence of such superposition 
rules is an interesting open problem. 

Proposition 2. Every system of SODEs (5) admitting a free superposition rule is a SODE 
Lie system. 

Proof. Suppose that system (5) admits a superposition rule of the special form 

X* = $^(xi,...,Xm;Ai,...,A2n), i = l,...,n. (10) 

In such a case, the general solution, x{t), of the system could be expressed as 

x'{t) = <^i,{xi{t),...,x„,{t);Xi,...,\2n), i = l,...,n. (11) 

If we set p{t) = {xi(t), . . . , Xm(t),xi{t), . . . , Xmit)), v^ = X*, and f^ = xjj with a = 1, . . . ,m, 
the derivative of the above expression with respect to t reads 

m n . „ j ^ 

v\t) = x^{t) = Y,T.U{t)-f{p{t))) , z = l,...,n. (12) 

a=lj=l ^ '^^« ^ 

Consequently, there exists a function 



^^{Xi,...,Xm,Vi,...,V„i) = 2_^l^ 



a=\ j=l ^ ^^« ^ 



vl — T , i = 1,... ,n, 



such that 

( x\t) = <^i{xi{t),...,Xmit);Xl,...,X2n), . ^ 

\ Z ^ 1 ... 71. 

[v'{t) = ^l{xi{t),...,Xrri{t),Vi{t),...,Vm{t);Xl,...,X2n), 

Therefore, system (7) admits a superposition rule and (5) becomes a SODE Lie system. D 

Apart from the SODE Lie system notion, there exists another method to study cer- 
tain systems of second-order differential equations which admit a regular Lagrangian, like 
Caldirola-Kanai oscillators or Milne-Pinney equations [10, 24]. Although this method can- 
not be used to study general second order systems, it provides us with some additional 
information that cannot be derived by means of SODE Lie systems when it applies, e.g. 
about the time-dependent constants of the motion of the system [24]. 

A possible generalisation of the concept of superposition rule appearing, for instance, in 
the theory of quasi-Lie schemes is the so-called time- dependent superposition rule [4]. This 
concept can be extended to the framework of systems of second-order differential equations 
as follows. 

Definition 3.3. We say that the map ^ : M x TM™" x M^*^ — ^ M" is a time- dependent 
superposition rule for the system of SODEs (5), if its general solution x{t) can be written 
in terms of each generic family xn-j{t), . . . ,X(jyi\{t) of particular solutions, their derivatives, 
a set of 2n constants, and the time as 

x{t) = '^{t,Xi(t),...,Xm(t),Xi(t),...,Xm(t);Xi,... ,A2n)- 

It is essential to analyse the existence of time-dependent superposition rules in order to 
understand the relevance of the practical results carried out throughout this work. As it 
shall be shown soon, many of the properties of these superpositions are a consequence of the 
features of time-dependent superposition rules for first-order systems. For instance, let us 
prove the following results concerning the existence of time-dependent superposition rules 
for systems of SODEs. 

Proposition 3. Every system of SODEs (5) admits a time- dependent superposition rule of 
the form ^ : M x M^n ^ ]^n_ 

Proof. Every system (5) is related to a first-order system (1) admitting a fiow g : (t; A) G 
M X R^" I—)- gt{X) G TM" which allows us to cast its general solution, ^{t), into the form 
^(t) = gt{k). Consequently, the general solution, x{t), of system (5) can be written as 
x{t) = IT o g^i^k). In other words, system (5) admits a time-dependent superposition rule 
depending just on 2n constants. D 

Despite the remarkable theoretical interest of the above result, it does not provide any 
additional method for the explicit derivation of solutions of systems of SODEs. Indeed, note 
that the derivation of the above time-dependent superposition rule amounts to working out 
the generalised flow for the flrst-order system (1). This involves solving the system for each 
initial condition. If this can explicitly be done, determining the above superposition rule 
becomes unnecessary; and, otherwise, the superposition rule, although interesting, cannot 
be provided. 



Apart from the time-dependent superposition rule guaranteed by Proposition 3, other 
instances can be ensured to exist. Nevertheless, their explicit determination uses to be as 
difficult as solving the initial system. Let us illustrate this statement more carefully. Con- 
sider a system of SODEs (5) related to the system of first-order differential equations (7), 
with general solution Px{t) = {x{t),v{t)). Let X be the time-dependent vector field associ- 
ated with this system, and Y any other time-dependent vector field on TM" ~ M^'". Their 
corresponding flows, g^ jQ^^ : M x M^" — )• TM", satisfy that Px{t) = gf- o[gY)~^{py{t)), where 
Py{t) is the general solution of the system describing the integral curves of Y. Therefore, 
(g^ o {g^)-'^)^Y = X (cf. [4]). In particular, if 1" = 0, then py{t) = (/)(Ai, . . . , A2n), where 
(f) : R^"' — ;• TM" is any diffeomorphism. Therefore, 

Px{t) = fff (Ai,...,A2„) =^ x{t) = vro^rf (Ai,...,A2n)- 

In other words, system (5) admits a time-dependent superposition rule depending just on 
a set of 2n constants. In a similar way, if we assume Y to be any autonomous, non- 
null, vector field, the system describing its integral curves is a Lie system ([1^,1^] = 0) 
and the straightforward application of the method developed in [4] shows that it admits a 
superposition rule depending on one particular solution Py^{t) = {yi{t),yi{t)), i.e. Py{t) = 
^{py^{t)\ Ai, . . . , A2n)- Consequently, the general solution, x{t), of the system (5) associated 
with X now reads 

X{t) =7Togfo (/)-! o $(gf o {gf)-\px^{t))- Ai, . . . , Asn), 

where Pxi{t) = {xi{t),xi{t)) is a particular integral curve of X. That is, system (5) admits 
a time-dependent superposition rule depending on one particular solution. 

Note that the determination of the above superposition rules relies, among other things, 
on the determination of the flow of the initial time-dependent vector field. Therefore, 
obtaining the superposition rule for (5) is as difficult as solving the system (5). Other 
similar constructions can be found. Nevertheless, most of them share the same drawbacks. 

The above remarks make evident that a key point in the study of time-dependent super- 
position rules for systems of SODEs is the development of procedures to explicitly determine 
them. As it occurred in the explicit determination of time-dependent superposition rules 
for systems of first-order differential equations, quasi-Lie schemes play an important role 
in the description of superpositions for systems of SODEs. Indeed, the same commentaries 
pointed out at the end of Section 2 can also be claimed here. 

On one hand, given a system of SODEs (5), quasi-Lie schemes provide a powerful tool of 
large applicability to derive time-dependent superposition rules for its associated first-order 
system (7). From here, it is immediate to obtain a time-dependent superposition rule for 
(5). 

On the other hand, quasi-Lie schemes naturally provide a family of first-order systems, 
including (7), whose elements admit the same obtained time-dependent superposition rule. 
It is easy to prove that all those systems of SODEs whose associated first-order systems are 
members of this family share the same time-dependent superposition rule. In Section 5 this 
fact will be clarified and illustrated through the study of second-order Riccati equations. 
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4 A new superposition rule for a family of SODE Lie systems 

In this Section we derive a superposition rule for a family of second-order differential equa- 
tions including, as particular instances, some Painleve-Ince equations [22]. In the process 
of searching for such a superposition rule, we find a family of Lie systems which will be used 
posteriorly to describe time-dependent and time-independent superposition rules for other 
second-order differential equations studied in Physics and Mathematics. 
Consider the family of differential equations 

x + 3x± + x^ = f{t), (13) 

with f{t) being any time-dependent function. The interest in these equations is motivated 
by their frequent appearance in the Physics and Mathematics literature [14, 15, 31]. The 
properties of these equations have been deeply analysed since their first analysis by Vessiot 
and Wallenberg [42, 44] as a particular case of second-order Riccati equations. For instance, 
these equations appear in [27] in the study of Riccati chain. In that work it is stated that 
the equations of the Riccati chain can be used to derive solutions for certain PDEs. In 
addition, equation (13) also appears in the book by Davis [18], and the particular case with 
f{t) = has recently been treated through geometric methods in [7, 14]. 

The results described in previous sections can be used to study differential equations 
(13). Let us first show that the above differential equations are SODE Lie systems and, 
in view of Proposition 1, they admit a superposition rule that is derived. According to 
definition 5, equation (13) is a SODE Lie system if and only if the system 

{X = V, 
(14) 
v = -^xv-x^ + f{t), 

determining the integral curves of the time-dependent vector field of the form 

Xpi{t,x,v) = Xi{x,v) + f{t)X2(x,v), (15) 

with 

r\ fi f^ 

Xi=v—-{3xv + x^) — , X2 = —, 
ox ov ov 

is a Lie system. 

In view of the decomposition (15), all equations (13) are SODE Lie systems if the vector 

fields Xi and X2 are included in a finite-dimensional real Lie algebra of vector fields V . This 

happens if and only if Xi, X2 and all their successive Lie brackets, i.e. the vector fields of 

the form 

[Xi,X2], [Xi, [Xi, X2]], [X2, [Xi, X2]], [Xi, [Xi, [Xi,X2]]], . . . (16) 

span a finite-dimensional Lie algebra. Consider the family of vector fields on TM given by 
ox ov ov 



A3 = --- + 3a;--, A4 = X- 2x --, 

Ox ov ox ov 

X^^{v + 2x')^ - x{v + 3x2)|-, Xe = 2x{v + x')^ + 2{v^ - x^)-^, 

d d d d 

X7 ^ x^-, Xs = 2x^— +'iv^—, 

ox ov ox ov 
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(17) 



where X3 = [Xi.Xa], -8X4 = [X^.X^^], X5 = [X^^X^], Xg = [X^^X^], X, = [X2,X5], Xg = 
[X2,X6]. Then, the vector fields Xi, . . . ,Xs are hnearly independent over M. Additionahy, 
in view of the previous commutation relations and 



[Xi,Xe]^0, [Xi,X7] = ^Xs, [Xi,Xs]^-2Xi, [X2,X3]=0, 

[X2,Xi]^0, [X2,Xr]=0, [X2,Xs]=4X2, [X^^X^^-Xj, 

[X^jXq] ^ —^Xs, [XsjXq] = —2Xi, [XsjXy] = — 2X2, [X3,Xg]^2X3, 

[X4,X5] = —Xi, [X4,Xq] ~ 0, [X4,X7] = X3, [X4,X8] = 0, 

[X5,X6]=0, [X5,Xr] = -3X4, [X5,Xs] = -2X5, [Xe,X,] ^ ^2X5, 

[Xe,Xs] = -4Xe, [Xr,Xs]=2Xy, 



(18) 



it follows that the vector fields Xi , . . . , Xg span an eight-dimensional Lie algebra of vector 
fields V containing Xi and X2. Therefore, equation (13) is a SODE Lie system. Moreover, 
the elements of the following family of traceless real 3x3 matrices 




satisfy the same commutation relations as the vector fields Xi, . . . ,Xs, i.e. the linear map 
p : s[(3,M) — >• V, such that p{Ma) = Xa, with a = 1,...,8, is a Lie algebra isomor- 
phism. Consequently, the finite-dimensional Lie algebra of vector fields V is isomorphic 
to s[(3,M) and the systems of differential equations describing the integral curves for the 
time-dependent vector fields 

8 
X{t,X,v)=^ba{t)Xa{x,v), (19) 

a=l 

are Lie systems related to a Vessiot-Guldberg Lie algebra isomorphic to s[(3, M). 

Recall that, in view of Proposition 1, a superposition rule for (13) can be obtained by 
means of a superposition rule for the Lie system (14). As we already stated in Section 2, a 
superposition rule for a certain Lie system describing the integral curves of a time-dependent 
vector field X admitting a decomposition of the form (19) can be obtained through the 
first-integrals for the distribution D spanned by the diagonal prolongations Xi . . . , Xg on 
a certain TM"^™"*"^) in such a way that their projections 7r*(Xa), with a = 1, ... ,8, are 
linearly independent on a dense open subset of TR"™. 

Consider the functions Fabc, with a, 6, c = 0, . . . , 4, given by 

Fabc = Va{Xc - Xb) + Vb{Xa - Xc) + Vc{xb - Xa) + {Xa - Xb){xb - Xc){Xc - Xa). (20) 

Such functions satisfy that Fabc = Fbca = F^ab = -Fbac = -Fcba = -Facb, for ah a,b,c = 
0, . . . , 4, and they are useful to determine the superposition rule we are looking for. If we 
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take now the diagonal prolongations to TM^ of the family of vector fields (17), we can get, 
e.g. by means of any symbolic manipulation program, that the vector fields n^Xa), with 
a = 1, . . . , 8, on TM^ are linearly independent at those points p = (xi, . . . , X4, vi, . . . , v^) S 
TM'* satisfying Fi23(p)-Fi24(p)-^i34(p)-^234(p) 7^ 0. Such a set is dense in TIR^ and hence 
the involutive distribution V spanned by Xi . . . , Xg on TM^ is eight-dimensional almost 
everywhere. Additionally, there exist, at least locally, two first-integrals for the vector fields 
of the distribution V which can be used to derive a superposition rule. 

As the vector fields Xi, X2 and their successive Lie brackets, i.e. the vector fields (16), 
span the distribution P, it can be proved that giving a first-integral F : TM^ — ?• M for the 
vector fields Xi and X2, i.e. XiF = X2F = 0, is equivalent to giving a first-integral for the 
distribution T). 

These integrals can be obtained by applying the method of characteristics to the vector 
fields Xi and X2. The calculation of these first-integrals is long and it is detailed in the 
Appendix. As a result, we get that two first-integrals for Xi and X2 and hence for all the 
vector fields of the distribution T) are 

Ai[p) = - — — - — — and A2[p) = - — ^-—=, — ^ 

F42l{P)F3lo{p) F42i{P)F4sq{p) 

where now p = (xq, xi, . . . ,Xi,vo,vi, . . . , Vi) G TM^. We can obtain vq from the expression 
of Ai as 

_ (r^l(x2 - Xq) + t;2(xo - Xi) + {xi - 3;o)(xo - X2)(X2 - Xi))F43i 
{X2 - Xi)F43i + (Xi - X3)F42lAi 
(r;3(xi - Xq) + 1)1 (xq - X3) + [xq - Xi)(xi - X3)(X3 - Xo))F42lAi 
(X2 - Xi)F43i + (Xi - X3)F42lAi 

and if we substitute this value in the expression for A2, we obtain the expression 

_ 3^2-^431 - G'3i24A2 — G'2134Ai + a;3f42lAiA2 , , 

^431 + (-^124 — -F324)Ai + (i^4i2 — -F312)A2 + AiA2i^421 

where 

Gabcd = Xa{{Vd - Vc)xb + {vb - Vd)Xc + (xfc - Xc)xbXc + (Xc - Xb)XaXd) + 

Xd{{Vc - Va)Xb + {Va - Vb)Xc + (Xc - Xb)XbXc + (Xf, - Xc)XaXd). 

Note that the first-integrals Ai and A2 satisfy that 

Aj(xo(t),...,X4(t),t;o(t),... ,V4:{t)) = cte. j = 1,2, 

where {xa{t),Va{t)), with a = 0, ... ,4, are any family of five particular solutions of system 
(14). Consequently, if we fix xo(t) = x(f), 'L'o(t) = xo(*) and take into account expression 
(22), we see that the general solution of system (14), namely, x(t), can be written in terms 
of each generic family of four particular solutions Xa{t), with a = 1, . . . , 4, their derivatives 
in terms of t, and two real constants Ai, A2 as follows 

, . ^ X2{t)Fi3lit) - G3124(0A2 - G2Mt)Xl + X3{t)F^2lit)XiX2 

^^' F43l(0 + (^124(t)-F324(t))Al + (F412(t)-F312(t))A2+AlA2F42l(i)' ^ ' 
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where the time-dependent functions Fabcif) and Gabcdif) are obtained evaluating the func- 
tions Fahc and Gabcd on the curves (xi(t), . . . , a;4(t),xi(t), . . . ^x^{t)). 

In order to ihustrate the previous result by means of a simple example, let us consider, 
for instance, equation (13) with f{t) = 0. By direct inspection, we can get the set of 
particular solutions 

2 2t 1 + 2t 

rEi(t)=0, X2{t) = -, X3it) = j^, X4(i) = y-p^. (24) 

These particular solutions can be used to determine the functions appearing in the expression 
(22), i.e. G2134, -F431, etc. More specifically, we get, in view of the particular solutions (24), 
that 



(^3124 (i) = /,9 , i\/, , i\ i -^431(0 = /,o , -|W, , -|^ ^ G2134(t) 



and 



(t2 + l)(i + l)' ^^^^^ (i2 + l)(i + l)' — V^^ (i2 + l)(i + l) 



2 2 2 

^124 (*) = ,9 , .3' ^324 (t) = ,9 , .3 , .4 , ..^' -^312 (i) 



^2+^3' -•^^^VV ^2+^3 + ^4 + ^5' -^^^VV 2t + t3' 

Finally, making use of expression (22) and the above functions, we get that the general 
solution for equation (13) is 

.it)- (i + 2a,)(-i + A2) 



i(-l + A2) + t2Ai(-l + A2) + (-1 + Ai)A2 ■ 

In view of Proposition 1, expression (23) is not only a superposition rule for the equation 
(13), but also for any SODE Lie system (5) satisfying that its corresponding system (7) is 
related to a time-dependent vector field that put into a form similar to (19), i.e. it takes 
values in V. Many instances of the family of Lie systems (19) are associated with interest- 
ing SODE Lie systems with applications to Physics or related to interesting mathematical 
problems. In all these cases, the theory of Lie systems can be applied to investigate these 
second-order differential equations, recover some of their known properties, and, possibly, 
provide new results. Let us illustrate this assertion by means of a few examples. 

Another equation appearing in the Physics literature [15, 16, 40] which can be analysed 
by means of our methods is 

X + 3xx + x^ + Xix = 0, (25) 

which is a special kind of Lienard equation x + f{x)x + g{x) = 0, with f{x) = 3x and 
g{x) = x^ + Xix. The above equation can also be related to a generalised form of an Emden 
equation occurring in the thermodynamical study of equilibrium configurations of spherical 
clouds of gas acting under the mutual attraction of their molecules [19]. 

As in the study of equations (13), by considering the new variable v = x, equation (25) 
becomes the system 

{X = V, 
■ \ 3 , (26) 

V = —6XV — X — \iX, 

describing the integral curves of the vector field X = Xi — Xi/2(Xi + X^) included in the 
family (19). Consequently, the expression (22) can be used to derive the solution of equation 
(25) in terms of a set of four particular solutions. 
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Finally, we can also treat by our methods the equation 

x + ?,xx + x-^ + f{t){x + x^) + g{t)x + h{t) = 0, (27) 

containing, as particular cases, all the previous examples [31]. The system of first-order 
differential equations associated with this equation reads 

{X = V, 
(28) 
i, = -3xv -x^ - f{t){v + x^) - g{t)x - h{t). 

Hence, this system describes the integral curves of the time-dependent vector field 

Xt = Xi- h{t)X2 - \f{t) {Xs - 2X4) - ^g{t) {X7 + Xs). 

Therefore, equation (27) is a SODE Lie system and the theory of Lie systems can be used 
to analyse its properties. In particular, the expression (22) provides us with the general 
solutions for these equations out of each generic set of four particular solutions. 

Some particular cases of system (27) were pointed out in [16, 31]. Additionally, the case 
with f{t) = 0, g{t) = uj'^{t) and h{t) = was studied in [15] and it is related to harmonic 
oscillators. The case with g{t) = and h{t) = appears in the catalogue of equations 
possessing the Painleve property [30]. Finally, our result generalises Vessiot's result [43] 
describing the existence of an expression determining the general solution of system like 
(27) (but with constant coefficients) in terms of four of their particular solutions, their 
derivatives and two constants. 

5 Quasi-Lie schemes and second-order Riccati equation 

In this Section we derive a time-dependent superposition rule for the second-order Riccati 
equation [14] 

X + (6o(t) + hi{t)x) x + ao(t) + ai{t)x + a2{t)x^ + a:i{t)x^ = 0, (29) 

with a3(t) > 0, 03(0) = 1, 6o(i) = -^=y - ^ and 6i(t) = 3^^^^, by means of the 

theory of quasi-Lie schemes. 

We first introduce a new variable v = x and transform equation (29) into the system of 
first-order differential equations 



\x = v, 

\v = — {bo{t) + bi{t)x) V — ao{t) — ai{t)x — a2{t)x — 03(^)2; . 

Consider the following set of vector fields 

d_ 

dv' 
d_ 

dx'' 
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(30) 



ox 


ov 


yz = xv , 

ov 


n 


OV 


ov 


Y,-x^^ 

^7 — X _ , 
OV 


Ys 



spanning a linear space of vector fields V = (Yi, . . . ,!§)• The system (30) describes the 
integral curves of the time-dependent vector field 

Yt = Yi- bo{t)Y2 - h{t)Y3 - ao{t)Yi - ai(i)l5 - a2(t)>6 - a^i^, 

and, therefore, as Yt G V, for every t G R, we get that Y G y(]R). Let us denote adz{X) = 
[Z,X], for every vector fields X,Z. Hence, as 



ad?-3(y6) = ady3 o . . . o ady3(y6) = (-^)"+']^, 



d_ 



each Lie algebra of vector fields V containing the vector fields ^3 and Yq must include 
the above infinite family of linearly independent vector fields over M. Consequently, as Y3 
and Yq are contained in V, there exists no finite-dimensional Lie algebra of vector fields V 
including V and, therefore, equation (30) is not a Lie system. 

However, we can deal with such a system as a quasi-Lie system with respect to a quasi- 
Lie scheme [4, Theorem 4] and therefore this could be used to build up a time-dependent 
superposition rule for equation (29). In fact, define the linear space W = (125^8)) which 
is a two-dimensional Abelian Lie algebra of vector fields and, in view of the commutation 
relations 



[^2,1^1] 


= Yu 


[12,13] 


= 0, 


[1^2,14] 


= -n. 


[Y2,Y,] 


= -n 


[^2,^6] 


= -Yq, 


[12,17] 


= -l7. 


[1^8,11] 


= -1^1, 


[l8,l3] 


= 1^3, 


[^8,n] 


= 0, 


[1^8,15] 


= 15, 


[l8,l6] 


= 2Yg, 


[l8,l7] 


= 3^7, 



we get that [VF, y] C V. Hence, the pair (VF, y) becomes a quasi-Lie scheme S{W,V). The 
theory of quasi-Lie schemes shows that the time-dependent vector field Y can be transformed 
through any element g of the group of transformations of the scheme G{W) into a new time- 
dependent vector field g^Y G ^(M), see [4, Proposition 1]. 

The time-dependent transformations associated with elements of ^(VF) are 

fx(i) = 7(t)x'(i), 

\v{t) = mv'{t), ^ ^ 

with 7(f) > 0, /3(t) > and 7(0) = /3(0) = 1. The above family of time-dependent changes 
of variables transforms system (30) into 

dt j{tf j{tf' 

dv' ao{t) .,^/^, .,w^aiW\ / Q2(t)7'(0 ,2 ,oo^ 

— 7(i) bi{t)v + — — r- X — — — X (32) 



dt (3{t) '^ ' y^ ' Pit) J p(t) 

_ a3(t)7'(t) ,3 _ &o(t)/3(t) + /3(t) 

m m 

In order to relate this system of first-order differential equations to a SODE, we have to 
choose 7(t) = K for a certain real constant K. As 7(0) = 1, then we choose 7(t) = 1 and 
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the previous system becomes 



dt 
It 



ao{t) 

m 



M.v.^).' 



m 



(33) 



a^it) ,3 6o(t)/3(t)+/3(t) , 
-X 7rr\ V . 



m 



m 



Let us try to relate this system to one of the Lie systems of the family (19), e.g. the one 
describing the integral curves of a time-dependent vector field Xt = fi{t)Xi + f2{t)X2 + 
{h{t)/2){Xs + Xr) + {U{t)/A){Xs - 2X4). If we fix /3(t) = y/^^ in (33), we obtain 

dx' 




ao{t) 



y^a^{3v'x' + X- 



/3\ 



ai(i) ^^, 



xAsW xAsW 



"^^'^ (.' + x'2), ^''^ 



and, consequently, the above system is a Lie system associated with the time-dependent 
vector field 



Xt = VasitiXi 






(X3 + X-r 



«2(i) 



4Va3(t) 



{Xs — 2X4 



As the above time-dependent vector field belongs to the family (19), the first component of 
the general solution (x' (t) , v' (t)) of system (34) can be described by means of the expression 
(23) in terms of four generic particular solutions {x'^{t) , Va{t)) . Moreover, as j{t) = 1 and 
I3{t) = Y^a3(t), then x{t) = x'{t), Xa{t) = x^(t) and v'^{t) = a^{t)-'^/'^dxa/dt{t). Using 
these relations in expression (23), we get a time-dependent superposition rule for second- 
order Riccati equations. More specifically, the general solution x = x{t) for any instance 
of second-order Riccati equation of the form (29) can be written in terms of a set of four 
particular solutions xi = xi(t), X2 = X2{t), X3 = X3{t) and X4 = X4(t) and its derivatives 
xi = xi{t), X2 = X2{t), X3 = x^it) and ±4 = X4(t) as 



2^2^431 — G2134A1 — G3124A2 — a;3^412AlA2 
^431 + (-^124 — ^324)Al + (F412 ~ ^312)A2 + A1A2F42I 



(35) 



where 

^abcd = Xa{a 



-1/2/ 



-1/2, 



3 {t){Xd - Xc)xb + 03 '''''{t){xb - Xd)Xc+ 

— 1/2 
{Xb - Xc)XbXc + {Xc - Xb)XaXd) + Xd{a^ {t){Xc - Xa)xb+ 

— 1/2 
03 {t){Xa - Xb)Xc + {Xc- Xb)xbXc+ {Xb - XcjXaXd), 



and 

^ abc 



-1/2 



Xa\Xc 



Xb) +a^ Xb{Xa 



Xc) + 



+ a. 



-1/2 
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with a,b,c,d = 1, . . . ,4. 

Finally, it is just worth noting that, as pointed out in Section 3, time-dependent superpo- 
sition rules for a system of SODEs appear naturally related to families of systems admitting 
the same superposition rule. Indeed, note that if we had initiated this section by analysing 
system (29) with as{t) = 1, we could have proceeded exactly in the same way as before. 
Nevertheless, when reaching the system (34) with a^i^t) = 1, we could have noticed that all 
systems with a more general 03 (t) are quasi-Lie systems admitting the same superposition 
rule as our particular instance. Consequently, this would have shown the existence of a more 
general system of SODEs, namely, (29), admitting the same superposition rule. 

6 Appendix 

We have relegated to this appendix various calculations that, although necessary to obtain 
certain previously stated results, i.e. the common first-integrals for vector fields of P, do not 
deserve to be detailed in the main body of the article as they do not provide any relevant 
knowledge. 

Recall that a function F : TM^ — )■ M is a common first-integral for every vector field in 
V if and only if XiF = X2F = 0, where Xi,X2 are the diagonal prolongations to TM^ of 
the vector fields Xi,X2. Therefore, such a function F must be a solution of the equation 

a=0 

written using the coordinate system {xq, . . . ,X4,vo, . . . ,^4}. This equation can be solved 
using the method of characteristics. Such a method explains that the solutions of the above 
equation are constant along the solutions, the so-called characteristics curves, of the system 

dxQ = . . . = dx4 = 0, dvQ = . . . = dvi. 

So, the characteristics of equation X2F = are curves {xq, . . . ,X4^,vq{s), . . . ,V4,{s)) such 
that ^0 = vq{s) - Vi{s), ii = vi{s) - va{s), S,2 = V2is) - va^s) and ^3 = ^3(5) - v^is) 
for certain real constants ^O) • • • ^^3- Thus, there exists a function i^2 : I^^ — ^ I^ such that 
F{xo, . . . , X4, vq, . . . , V4) = F2{xo, . . . , Xi, ^0, . . . , ^3). In other words, function F depends 
only on the variables xq, . . . , X4, ^Qj ■ ■ ■ ; Cs- 

Consider the coordinate system {xq, . . . , X4, ^O) • • • ) Cs) ^4} on TM^. In terms of the new 
coordinate system, the vector field Xi reads 
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^1 = E [^-^^ - (3^-^^- - ^4 + 4) ^) - (4 + 2>xaVa)^+ 



As we assumed that F is a first-integral of the distribution P, it is a first-integral for the 
vector fields Xi and X2. Taking into account that F is a solution of the equation X2F = 
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and, therefore, it depends only on the variables xq, • • • , x^, Co) ■ ■ ■ ; Cs; the equation XiF = 
yields 



t 0»f - (-"«» -t . 4) f ) . n (e f - 3 E(. - .,f ) 



0, 



and in view of the dependence of the function F2 , we get that there exist two vector fields 



a=0 



0=0 



a=0 



such that XiF = Z1F2 + V4Z2F2 = 0. Therefore, we have that Z1F2 = and Z2F2 = 0. 
In consequence, the function F is a first-integral of the distribution T> if and only if it is a 
first-integral for the vector fields Zi and Z2 depending on the variables xq, . . . ,Xi,S,o, ■ ■ ■ ,^3- 
The first-integrals of the vector field Z2 depending just on the above variables can be 
determined by means of the characteristic curves given by the following system 



d{Xa - X4) = 0, 



dxA 



Z{Xa - X^) 

Hence, such first-integrals depend on the functions 



a = 0, 



,3. 



'Ha, — Xa X4., 
(t>a = "iflaXA + ia-, 



a = 0, ... ,3, 



Consequently, given a first-integral F of the distribution P, there exists a function F3 : 
M^ — ;• M such that F{xq, . . . , X4, vq, . . . , V4) = ^3(770, . . . , ??3, (^O) • • • ; 4*3)^ i-6- the function F 
actually only depends on the variables %, ...,%, 0O) •••) '/'a- 

Taking into account the dependence of the function F in the above variables, the equation 
ZiF = reads in the coordinate system {rjo, ...,7/3, cpQ, . . . , cp^, X4, ^4}, 



E 

a=0 



OVa 0<Pa 



0. 



Collecting terms with different powers of 3:4, we obtain that Z2F = ^qFj, + 3:4^1^3 
3x1^2^3 = 0, with 



f^n 



^^„ . "Wa - ^''^"^^ + ^'^ a^) ' ""'= S {-''^Wa + '^'a^) ' 



E 



02 



a=0 



d 



As F3 does not depend on 2:4 in the set of coordinates we have chosen, then f^a-^s = ^aF = 0, 
for a = 0, 1, 2. The method of characteristics for the equation i^iF^ = implies that there 
exists a function F4 : R'^ — )• M such that F{xq, . . . , X4, vq, . . . , V4) = F^^Sq, . . . , (^3, Li, L2, L3) 
with 

' 5a = (pa + vl, a = 0, ... ,3, 



r _ '^'^ 

-f^a — — ) 



a = 1,2,3. 
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Now, using the coordinate system {60, . . . , ^3, Li, L2, Ls,r]Q, X4, v^}, we get 



n..,.n.,,.,,(2^^,±yJiyo. 



and, repeating the previous procedure, we see that there exists a function i<5 : M^ — )• M 
such that F(xo,... ,X4^,vo,... ,V4) = F5(Li, L2, L3, Ai, A2, A3), where A^ = La6o - 5a and 
a = 1,2,3. 

As we have shown that finding a first-integral F for the vector fields of the distribution 
T> reduces to looking for a first-integral F5 of the vector field f^o depending on the vari- 
ables Li, L2, -^3, Ai, A2, A3, we still have to analyse the condition VLqF^ = to determine 
completely the form of the first-integrals for the distribution T). 

By choosing the coordinate system {Li,L2, L3, Ai, A2, A3, (5o,r/o,X4,f4}, the equation 
O0-F5 = reads 



.^t(<^"-^"'t-^-"S)-(t""g-»-"S)-' 



3 
^0' 



a=l ^ \a=l 

and considering the vector fields 

a=l a=l a=l 

and the form of the function F5, the equation Oq-Fs = implies that Hii<5 = H2-F5 = H3F5 = 
0. 

If we apply the method of characteristics to the equation H3F5 = 0, it yields that there 
exists a function Fg : M^ — )■ M such that -F(xo, • • • , Xi,vo, . . . , v^) = Fq{Li, L2, L3, 712, 713), 
with TT2 = A2A^ and tt^ = A3Aj~ . Moreover, as F^ also satisfies the equation H2-F5 = 
^2-^6 = 0, we obtain that Fq only depends on the variables tt2, its, T2 = t^2Li — L2 and 
r3 = TT^Li — L^, i.e. there exists a function F7 : R^ — ;• R such that F{xq, . . . , X4, vq, . . . , v^) = 
Fe{Li,L2,L3, 712,113) = Fj{iT2,7r3,T2,T3). 

Finally, the conditions HiF = 'EiFj = imply that H1F7 = T2-F7 + L1T1F7 = 0, where 

2.aF7 , , 2^5^7 ^ dFr ^ OFj 



and 



^ ^ 2^57r2 ^ ^^57r3 5r2 ar3 



T2F7 = (r2 + r2)^ + (rs + r3)^ + r2vr2— + r3vr3— = o. 



The first equation implies that there exists a function Fg : R"^ — ;• R which satisfies that 
F{xo,...,Xi,vo,...,Vi) = F7{7r2,-K3,T2,T3) = F8(^o,*i,*2), where *o = (7r2 - l)r2"\ 
^1 = ("^3 — l)r3~^ and ^2 = T^2'^3~^{^ — 7r3)(l — 1^2)'^- So, the equation T2-F7 = T2-F8 = 
can be cast into the form 

dFj dFr ^1 - 'I'n dFr 

^ ^'d^O ^ ^9*1 ' *o*l 9*2 
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and we obtain that F is an arbitrary function of the functions Ai = j—^ and A2 = ^ ^ . 
In particular, undoing ah the changes of variables performed along this Section, we get that 

Klip) = - — --- — -- and A2 p = - — ^^-^^ — ^^ 
FA2i[p)Fua{p) F42l[p)F4so{p) 

are two independent first-integrals for the vector fields of the distribution T>. 

7 Conclusions and outlook 

We have defined and analysed the concepts of superposition rule, time-dependent superpo- 
sition rule, and free superposition rule for systems of SODEs. Several results concerning the 
existence of such superposition rules have been proved. Posteriorly, our theoretical achieve- 
ments have been illustrated by means of the study of a number of SODEs appearing in the 
Physics and Mathematics literature. Several new SODE Lie systems have been described 
and a common superposition rule for all of them has been derived. Such a superposition 
rule has posteriorly been used, with the aid of the theory of quasi-Lie systems, to get a 
time-dependent superposition rule for second-order Riccati equations. 

In the future, we expect to continue the analysis the properties of superposition rules 
for systems of SODEs as well as to investigate the generalization of the techniques depicted 
throughout this work to systems of higher-order differential equations. As an application, we 
hope to apply the results obtained in order to study new systems of second- and higher-order 
differential equations. In particular, it is specially interesting analysing the higher members 
of the Riccati hierarchy in other to develop new methods to determine solutions for those 
PDEs whose Backlund transformations are described by members of such a hierarchy. 
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